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Emergence of superconductivity in a two-band system coupled with breathing and Jahn-Teller phonons 
is discussed in a weak-coupling limit. With the use of a standard quantum mechanical procedure, the 
phonon-mediated attraction is derived. From the analysis of the model including such attraction, a BCS- 
like formula for a superconducting transition temperature T c is obtained. When only the breathing phonon 
is considered, T c is the same as that of the one-band model. On the other hand, when Jahn-Teller phonons 
are active, T c is significantly enhanced by the interband attraction even within the weak-coupling limit. 
Relevance of the present result to actual materials such as iron pnictides is briefly commented. 
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After the Bardeen-Cooper-Schrieffer (BCS) theory for su- 
perconductivity, 1 * it has been pointed out that anisotropic 
Cooper pairs can be formed by the Friedel oscillation in elec- 
tron systems, mainly from an academic viewpoint. 2 * Nowa- 
days it is widely recognized that anisotropic superconductiv- 
ity originating from strong electron correlation confirms an 
important route to achieve high superconducting transition 
temperature T c . In fact, due to successive discoveries of su- 
perconductivity in strongly correlated electron systems such 
as molecular conductors, transition metal oxides, and heavy- 
fermion compounds, it has been one of central issues in the 
research field of condensed matter physics to elucidate the 
mechanism of anisotropic superconductivity with relatively 
high T c . Among them, concerning a superconducting mate- 
rial group characterized by singlet Cooper pair, a key concept 
of d-wave superconductivity mediated by antiferromagnetic 
spin fluctuations has been believed to be established. 3 ~ 8) 

In addition to the concept of anisotropic superconductiv- 
ity mediated by magnetic fluctuations, another important in- 
gredient is multiband effect, since superconductivity has been 
found in electron systems with multiband such as S^RuCi, 9 ' 
MgB2, 10) and iron pnictides. 11 ' 12) Since multi-sheets of Fermi 
surfaces have been usually observed in heavy fermion com- 
pounds, the occurrence of superconductivity in such materials 
should be also related to multiband nature. Just after the BCS 
theory, multiband effect on T c has been discussed in a simple 
two-band electron model. 13) In fact, in recent years, multiband 
superconductivity bas been actively investigated from various 
viewpoints. 14 ~ 22 > In particular, here we mention s±-wave su- 
perconductivity proposed for iron pnictides. 23 ~ 26) 

In general, in electron systems with degenerate orbitals, 
Jahn-Teller phonons should play an important role, since 
Jahn-Teller distortions are known to lift the degeneracy in 
electron orbitals. In fullerene superconductors, 27 ~ 30) s-wave 
pair formation due to Jahn-Teller phonons has been discussed. 
A possibility of superconductivity due to geometric phase in 
Jahn-Teller crystals has been proposed. 31,32) However, the at- 
tractive interaction mediated by Jahn-Teller phonons has not 
been analyzed satisfactorily even in a weak-coupling limit, 
probably because of tedious calculations to derive such effec- 
tive attractions in multiband systems. 

In this paper, the effective interaction is derived in the two- 
orbital electron system which is coupled with breathing and 



Jahn-Teller phonons with the use of a standard quantum me- 
chanical techniques. By applying a weak-coupling approxi- 
mation to the effective model, we discuss the enhancement 
of T c due to Jahn-Teller phonons in the two-band electron 
system. When we include only the breathing phonon, T c is 
just the same as that of the famous BCS formula for the 
one-band model. On the other hand, when we consider Jahn- 
Teller phonons, we find the significant enhancement of T c by 
the interband attraction. These weak-coupling solutions are 
checked from the numerical estimation of pair susceptibility 
in the effective 4-site model. It is concluded that T c of the 
multiband systems coupled with phonons is enhanced by the 
number of relevant phonon modes. 

Now we consider a two-orbital electron system coupled 
with breathing and Jahn-Teller phonons, given by 



H = 



+a jt J. 



, gtQupu 



\Pl/{2M l ) + k t Ql/2], 



(1) 



where di 1<y is an annihilation operator for an electron with 
spin a in the orbital 7 (=a and b) at site i, t% , denotes the 
electron hopping between adjacent 7- and 7' -orbitals in near- 
est neighbor sites connected by a vector a, gi indicates the 
coupling constant between electrons and local distortion spec- 
ified by I, Qu indicates breathing distortion, Q 2 i and Q-^i de- 
note Jahn-Teller distortions, p u = d\ d iac , + d\ ha d iha , p 2 i= 



Aaa^ibo + d\ ba d iaa , p 3i - d\ a(J d iaa - d\ b<J ^ 

canonical momentum of Qn, Mi is corresponding reduced 
mass, and kg denotes the spring constant. 

By following the standard procedure of quantization of 
phonons, we introduce the phonon operator an, defined 
through Qn= (a 



d] h „diba, Pu denotes 



at •) / \/2u)tMg., where we is the phonon 



"IV 

energy given by wi=\fktjMg. After performing the Fourier 
transform, we obtain H in the form of H=Hq+Hi, where Ho 
denotes the sum of electron and phonon energy, given by 

Ha = E kT<^ kTIJ d hTa + ^ ue(a\ q ae q + 1/2). (2) 
Here the electron energy is given by 



Ekr — £kaa + Efefab ± sj (Cfcaa ~ £kbb) 2 + ^kab (3) 
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where r=l and 2 correspond to + and — signs, respectively, 
t a 



and £fc 77 '= J2 a e lfc a £" 7 " The electron-phonon coupling term 



Hi is given by 

Hi= > y/aiui 



- E 

e,kq,TT' 



,)C(fc,9)^(fc,?), (4) 



where ae=gj/(2M e ivf), p TT ,(k,q)= ^ d k+qTtT d kr >a, and 
the coefficient matrices T's are given by 

r (1 \k,q)-. 

f (2) (fe,9)=( , )■ (5) 

f (3) (fe,9) = 



Here A^_= u£, 



± fe,q _ u, k+q U k 

u k+ q u k> wnere is s iven b y 



+ 



± u 



A + 

k,q 
~ B k,q 


fl 5* ) 

A tq ) 


~ B k,q 
A k,q 


B L ) 


A k,q 
B k,q 


B L ) 

- A k,q J 


k+q U k 





l k + q U k ± 



1± 



£ka 



Skbb 



\/(Jka 



Skbb) 



^ b kab 



1/2 



(6) 



Since we assume degenerate the Jahn-Teller modes here, we 
set ix>i=(x>br and lj 2 =uj 3 =ujjt- Concerning coupling constants, 
we introduce ai=a^ r and a 2 =a^=ajT- 

Now we derive the effective Hamiltonian H c ff from H due 
to the elimination of phonon degrees of freedom by using 
a canonical transformation. Let us here consider a transfor- 
mation H e ff=e s He~ s , where the operator S is defined so as 
to satisfy the relation Hi=— [H , S]. Then, after some calcu- 
lations of operators, we obtain the effective Hamiltonian as 
H e f[=H + H int . The effective interaction between electrons 
mediated by phonons is given by 

H int = -[[H ,S],S)/2. (7) 

Here we obtain the effective model within the second order of 
electron-phonon coupling constant. 

In the present case, first we assume S in the form of 



X+ T ( ?\k,q)4_ Q }p TT ,(k,q), 



(8) 



where X^P is determined so as to satisfy Hi=—[H , S}. Af- 
ter lengthy calculations, we obtain 



X^\k,q) 



Ek+qr — Ekr' ± Wl 



i?>(M). 



(9) 



The effective interaction H- lnt is evaluated by [Hi , S] /2 from 
eq. (7). Then, we obtain 



^ ^l(k,q)T^,(k',- q ) 



H h 



E 2. 

i,TT' 'jjUju' k.k'.q 



(E k < T , - E k '_ qT ) 2 -u> e (io) 



x p T T<(k, q )Pw>(k',-q), 

where Ug=aguj£ and Ui denotes the attractive interaction 
due to £-mode phonons. Note here that U^U^^a^uj^ and 

^2=t^3=kjT=ajTWjT- 

Here we provide a comment on the Cooper pair in the sys- 




-71/2 



Fig. 1 . Cooper pairs in a two-dimensional electron system with a couple of 
Fermi surfaces. In the weak-coupling limit, pairs are formed by electrons 
on the same Fermi surface, since the pair between different Fermi surfaces 
should be composed of electrons with p 1 and — p 2 . 



tem with multi-Fermi surfaces. In the weak-coupling limit, we 
usually consider the pairing of electrons only in the vicinity 
of the Fermi energy Ep. Thus, we set E k < _ qT =E k < T ,mE-p. As 
shown in Fig. 1, in the weak-coupling limit, the Cooper pair is 
formed only by the electrons on the same Fermi surface, ex- 
cept for an unrealistic case in which a couple of Fermi-surface 
sheets are perfectly degenerate. Then, we consider the inter- 
action for Cooper pair with zero total momentum. Note that in 
the strong-coupling region, it is possible to consider the pairs 
of electrons far from the Fermi surface, leading to a chance of 
pair formation between different Fermi surfaces. 

After some algebraic calculations, we obtain the effective 
interaction in the form of 

ffint = - V ^(P'P') d l^ d - P ^ d -P'H d P'^- (11) 

P,p' 

The pair potential is given by 

Vu(p, P ') = V 22 (p,p') = U, T + U hl A;? p _ pl , 
Vi 2 (p,p') = V 21 (p,p') = U, T + U hr B p , 2 p _ pl . 



(12) 



Here we note the relation of A p ? p _ p ,+ B p , 2 p _ p ,=l. 

Next we solve the gap equation. The gap function A^p) 
is given by 



= ^y^{p,p'){ d -p'H d p'^), 

p' ,v 



(13) 



where (• • • ) denotes the average by using H in the mean-field 
approximation. Note again that it is enough to consider the 
pairs on the same Fermi surfaces in the weak-coupling limit. 
By assuming the Cooper pair with s-wave symmetry, we ob- 
tain the gap equation at T=T C as 



A M = log(1.13o; c /T c ) 2_, A^A„ 



(14) 



where uj c is an appropriate cut-off frequency and is the 
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Fig. 2. (Color online) Curves of (1) T c /u) c = 1.13e -1 /*br, (2) T c /w c = 
1.13e- 1 /< 2A JT), and (3) T c /ui c = 1.13e _1 /( 2AjT+A *"'). Here we set 
Ajx = Ab r = A for simplicity. 



non-dimensional coupling constant, given by 

All = A 22 = Ajt + Abr/3, 
Al2=A 2 l=Aj T + Abr(l-/3)- 



(15) 



Here Ajt=NqUjt, Abr = -^o£/br» Nq denotes the density of 
states at the Fermi level, (3=(A p ? ' p _ pl )FS> an d (• ■ -)fs de- 
notes the average over the Fermi surface. Note that we simply 
assume the same values of iVo for the different Fermi surfaces. 
By solving the gap equation eq. (14), we obtain 

T c = 1.13w c e" 1/(All+Al2) = 1.13w c e" 1/(2AjT+Abl) . (16) 

Note that another solution provides smaller T c even if it exists. 
Equation (16) tells us several interesting stories. First we con- 
sider a situation in which only the breathing mode is active. 
In this case, we immediately obtain the same formula of T c , 
r c =1.13w c e _1/AbI , as that of the BCS theory for a one-band 
system. Namely, even if the number of the band is increased, 
the magnitude of T c is not changed as long as the total electron 
density is coupled with the breathing phonons. In this situa- 
tion, there is no advantageous points of multi-band nature for 
the elevation of T c . 

On the other hand, in the multi-orbital system, there oc- 
curs a coupling with Jahn-Teller phonons so as to lift the de- 
generacy in electron systems. In such a case, the factor 2 ap- 
pears in front of the coupling constant in the T c formula. In 
other words, this factor 2 indicates the number of Jahn-Teller 
modes, not the number of electron bands. When we consider 
the coupling of degenerate electrons with both breathing and 
Jahn-Teller phonons, the factor 3 becomes effective if we sim- 
ply consider Abr=AjT- In Fig. 2, we show the curves of T c 
for the three cases of (1) breathing phonon, (2) Jahn-Teller 
phonons, and (3) both breathing and Jahn-Teller phonons. 
Here for simplicity, we set Abr=AjT=A. As easily understood, 
the change of the factor in the power is remarkable, even in 
the weak-coupling approximation. Note that the factor 2 or 3 
indicates the total number of phonon modes which are cou- 
pled with electron systems. 

In order to confirm the present result of eq. (16) obtained in 
the weak-coupling BCS approximation, we evaluate the pair 
susceptibility of the effective two-band model with the attrac- 
tive interaction induced by phonons. For the purpose, we re- 
sort to an unbiased technique such as exact diagonalization. 



The singlet pair correlation function is evaluated in a small- 
sized cluster and the effective coupling constant A c ff is de- 
duced from the singlet pair correlation. 
The effective model is given by 

H e s = ^2 E krdl Ta d kTa - I'Y] d iTf d iT fd iTl d iTi 

fe,T£T i,r 



3 ^2( d ln d in d in d m + h - c -), 



(17) 



where / and J denote on-site and pair-hopping attractive in- 
teractions corresponding to V\\ and V\%, respectively. In or- 
der to reproduce the situation in eq. (15), we set / and J as 
I = Ujt + Ubi/3 and J = U JT + Ub T (l —(f), respectively, 
by taking j3 as a parameter. 

The pair susceptibility matrix x Alt/ (m, n) is defined by 

X^(m,n)= dr(0 m , M (r)0t ), (18) 



-Ht 



, m m- 



where T is a temperature, (f> m ,n(T~)= e 
dicates the vector connecting possible two sites in the cluster, 
and (j>m.n is a singlet pair operator of the band /i, given by 



The coefficient . 



fim,v( d il4- d i+mtf — dirfdi+ m pi)\/2. (19) 



j is determined by the diagonalization of 
the susceptibility matrix and the pair susceptibility x^v is de- 
fined by its maximum eigenvalue. 

In order to extract information on the pairing interaction, 
we consider the pair susceptibility in a diagrammatic man- 
ner. When we define the non-interacting pair susceptibility as 
X^°\ X satisfies the relation of x = X^+ X^Kji m tne ladder 
approximation, where K denotes the effective pairing inter- 
action in the matrix form. Usually we obtain x from K and 
X^ \ but here we evaluate K as K=x^~ x — X" 1 * where x 
is numerically evaluated. Then, we can evaluate the effective 
non-dimensional coupling constant A e fr as 



Acff = (Kn + K 12 )/W, 



(20) 



where W is the bandwidth. 

In this paper, for the evaluation of x, we exploit an exact di- 
agonalization technique for the model in a 4-site cluster. We 
set£fci = — i(cos k x + cos k y ) and Ek2 — t(cos k x +cos k y ), 
leading to a couple of Fermi surfaces around V and M points 
in the thermodynamic limit. Here we consider the case of 
n=0.5, where n is the electron number per site and per orbital. 
In this case, X\\ -X22 -4/i and W=4t in the 4-site cluster. 

In Fig. 3(a), we show A e g vs. f3 by solid symbols for 
£/br=0.05i, O.lt, 0.2t, and O.Zt with n=0.5, [7 JT =0, and T=0. 
Note that each horizontal line denotes the pair susceptibility 
of a one-band model with nearest neighbor hopping t and 
on-site attraction Ubr for the same value of n in the 4-site 
cluster. For small values of Ubr/t, we find that A ff weakly 
depends on (3 and it agrees well with the result of the one- 
band model, suggesting that T c in the two-band system cou- 
pled with breathing phonons is just the same as that of the 
one-band model. Namely, in such a case, T c is not expected 
to increase even if we increase the number of electron bands. 

In Fig. 3(b), we show A e fi vs - U jt with a fixed value of 
[3=0.5 for three cases as (1) U^=U and L/jt=0, (2) t/br=0 and 
Ujt=U, and (3) Ub r =UjT=U . The effective coupling constant 
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Fig. 3. (Color online) (a) Effective coupling constant A e g vs. /3 for 
U hl /t=0.05, 0.1, 0.2, and 0.3 with n=0.5, U JT =0, and T=0. As for the 
horizontal lines, see the main text, (b) A e jj vs. U/t for (fbr> Ujt)=(U ,0), 
(Q,U) and (U,U) with /3=0.5. Three lines indicate A cff =0.022 X U/t, 
0.044 x C//i, and 0.066 x U/t. 



is found to be in proportion to U/t. Among the proportional 



coefficients, we confirm the relation of X^g = A^g /2= A^ /3, 
which is consistent with the weak-coupling result of eq. (16). 

Thus far we have analyzed the models in the weak-coupling 
region, but we are also interested in the strong-coupling be- 
havior. For instance, when we increase Ubr in Fig. 3(a), /in- 
dependence of A c ff becomes more significant and the devia- 
tion from the one-band result is large. When we increase the 
value of attractive interaction in Fig. 3(b), A e fj is no longer 
in proportion to U/t and the relation of A^L = A2jL /2= A2L /3 
does not hold. In order to discuss such strong-coupling ef- 
fects, we should analyze directly the original Hamiltonian 
eq. (1), not the effective model eq. (17), by applying the 
Migdal-Eliashberg theory. It is one of future problems. 

Here we provide a brief comment on the polaron mass en- 
hancement, which is one of strong coupling effects. For Hol- 
stein phonons, it has been well known that the polaron mass 
m* is increased as m* / m=e abI , where m is the bare electron 
mass, while for Jahn-Teller phonons, the mass enhancement 
has been found to be expressed by m*/m«e QjT / 2 / v /ajr for 
large ajx- Namely, for the same value of the coupling con- 
stant, the mass of the Jahn-Teller polaron is smaller than that 
of the Holstein polaron. This behavior seems to be related 
to the fact that the vertex corrections in an electron system 
coupled with Jahn-Teller phonons should be less effective in 
comparison with the case of Holstein phonons. 33 ' 34) The fact 
may be also relevant to the increase of superconducting T c in 
electron systems coupled with Jahn-Teller phonons. 

Note also that in the Migdal-Eliashberg theory, the effect 
of Coulomb interaction is included in the parameter fj,* as the 
reduced repulsion. The same story is expected to be applied 
to the present case, as long as we consider adiabatic phonons. 
However, the degree of the reduction of Coulomb repulsion 
may be different between intra- and inter-orbital interactions. 
If the on-site interaction is still negative while the pair hop- 
ping interaction becomes positive, we obtain the so-called s±- 
wave pairing, proposed for iron pnictides. 23 ~ 26) Thus, it may 
be possible to construct an alternative Jahn-Teller phononic 
scenario for superconductivity in iron pnictides. In particular, 



v(3) 



competition between Coulomb repulsion and phonon-induced 
attraction may be a key issue to understand the appearance of 
both nodal d-wave and nodeless s±-wave gaps in iron pnic- 
tides. This is an interesting problem in future. 

In summary, we have discussed the appearance of su- 
perconductivity in the two-band electron system coupled 
with breathing and Jahn-Teller phonons in the weak-coupling 
limit. It has been found that T c is increased with the increase 
of the number of relevant phonon modes. Namely, T c of the 
two-band system coupled with Jahn-Teller phonons becomes 
high in comparison with that of the one-band case, leading to 
a possibility of phonon-induced high-T c superconductivity. 

This work has been supported by a Grant-in-Aid for Sci- 
entific Research on Innovative Areas "Heavy Electrons" (No. 
20102008) of The Ministry of Education, Culture, Sports, Sci- 
ence, and Technology, Japan. 
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